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We investigate a model of closed {d — l)-dimensional soft-self-avoiding 
random surfaces on a d-dimensional cubic lattice. The energy of a surface 
configuration is given hy E = J{n2 + 4/0714), where 77-2 is the number of 
edges, where two plaquettes meet at a right angle and 774 is the num- 
ber of edges, where 4 plaquettes meet. This model can be represented 
as a Z2-spin system with ferromagnetic nearest-neighbour-, antiferromag- 
netic next-nearest-neighbour- and plaquette-interaction. It corresponds to 
a special case of a general class of spin systems introduced by Wegner and 
Savvidy. Since there is no term proportional to the surface area, the bare 
surface tension of the model vanishes, in contrast to the ordinary Ising 
model. By a suitable adaption of Peierls argument, we prove the existence 
of infinitely many ordered low temperature phases for the case A: = 0. A 
low temperature expansion of the free energy in 3 dimensions up to order 
x^^ {x = e~^'^) shows, that for k > only the ferromagnetic low tempera- 
ture phases remain stable. An analysis of low temperature expansions up 
to order for the magnetization, susceptibility and specific heat in 3 di- 
mensions yields critical exponents, which are in agreement with previous 
results. 



1 Introduction 



The so called gonihedric string was introduced by Savvidy et. al. [jlj-Q as a 
model for closed triangulated random surfaces. The action is given as 



s= Yl l^i-^jl^K 



simple edges 
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common edges 



The first term in (|I|) sums over all edges, where two triangles meet. The 
function 0{aij), which weights the edge lengths \Xi — Xj\, depends only on the 
angle aij between the neighbouring triangles such that 

1. 9{2n-a) = e{a) 

2. 9{n) = 

3. e{a) > 0. 

If more than two triangles meet at a given edge, contributions from all pairs 
of triangles arise in the second term of ([l|). Therefore if > 0, this term pe- 
nalizes self-intersections of the surface. Since 6{7i) = 0, the gonihedric action 
is subdivision invariant, i.e. geometrically nearby surface configurations have 
similar weights. The action measures essentially the linear size of the surfaces. 
Configurations with long spikes are therefore suppressed. These kind of con- 
figurations destroy the convergence of the partition function for triangulated 
random surfaces with area action p. However, for /c = and 9{a) = ||7r — q;|, 
the gonihedric action was shown to suffer from a similar disease 0. In this 
case fiat "pancake like" configurations dominate such that the grandcanonical 
partition function fails to converge. Nevertheless numerical simulations of the 
canonical ensemble |]8|-[10| show fiat surfaces. 



Another way to define discretized random surface theories is to consider pla- 
quette surfaces on a cubic lattice. In this approach, not just the surface, but 
also the embedding space is discretized. The gonihedric string can be formu- 
lated as a model for plaquette surfaces on a euclidean lattice as was shown 
by Wegner and Sawidy |[TTHT6|. If self-overlapping of the surface is excluded 
(i.e. each plaquette occurs either once or not at all in the surface), the model 
for closed {d — l)-dimensional plaquette surfaces can be written as a Z2 Ising 
model. The spins a = ±1 sit on a d-dimensional hypercubic lattice and the 
surface consists of all plaquettes of the dual lattice, which separates spins of 
opposite sign. The energy of a surface configuration on the lattice is now given 
as 

E=J2 ^cdgo = J{n2 + Akn^), (2) 

edges 



where n2 is the number of {d — 2)-dimensional edges, where two plaquettes 
meet perpendicular and n4 is the number of edges, where four plaquettes meet 
(see figure |l]). Expressed in terms of spin variables, the Hamiltonian reads 

k 1 — k 

H = -Jk{d -I) ^i^3 + ^ 9 XI ~ ^ H o-jo-jo-fca;. (3) 

<ii> [ijkl] 
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-^cdgc J -^cdgc -^cdgc 4:hJ 

Fig. 1. The total energy can be written as a sum over edge contributions. "Simple 
edges" , where two plaquettes meet perpendicular cost energy J. "Flat edges" , where 
two plaquettes of the same orientation meet, do not cost energy. If four plaquettes 
meet at a given edge ("double edge"), the corresponding energy contribution is AkJ, 
since according to the gonihedric string action four right angles are counted. 

It contains ferromagnetic nearest neighbour (< ij >), antiferromagnetic next 
nearest neighbour (<< ij ») and plaquette terms {[ijkl]), where the corre- 
sponding couplings are tuned in a special way. Similar surface models, where 
the energy contains an additional term proportional to the surface area, have 
been considered earlier |[T9|-23|, in particular in connection with amphiphilic 
systems. However, the special choice of couplings as given in (^, has not been 
studied explicitly in this context. It corresponds to the disorder line as calcu- 
lated in mean field approximation p2|. In two dimensions, the model can be 
mapped to the eight- vertex model |2^. The corresponding weights are how- 



ever different from the exact solvable case. Numerical simulations indicate, 
that the system defined in (H) does not undergo a phase transition for (i = 2 
but is in a disordered phase for all T > p5| , p!7| . In three dimensions the go- 
nihedric Ising model (0) has been studied by means of mean field I^G], Monte 
Carlo |[T7|j26| -|28| and cluster variation-Pade methods [^,0. At A; = the 
model seems to undergo a first order phase transition P^, whereas for large 
k the transition becomes second order [^,^,^. Since flat edges, i.e. edges 
where two plaquettes of the same orientation meet, do not cost energy, the 
ground state is highly degenerate. For > it consists of all configurations, 
where only flat domain walls are present, as long as they do not intersect. 
The ground state degeneracy of a finite system of size L"^ therefore increases 
like d2^. If /c = 0, the number of ground states is even higher, since ground 
state planes are now allowed to intersect without an additional energy cost. 
Thus the degeneracy increases like 2'^^ in this case. Moreover, the flipping of 
a whole {d — l)-dimensional spin layer does not change the energy of any spin 
configuration if = as can easily be seen from (^, i.e. apart from the global 
spin-flip symmetry the system possesses an additional layer-flip symmetry. 



The outline of the paper is as follows: In section 2 we show by a suitable 
extension of Peierls contour method, that the layer-flip symmetry for /c = is 
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spontaneously broken at low temperature. For each ground state we find an 
ordered low temperature phase. In section 3 we perform a low temperature ex- 
pansion around all possible ground states of the gonihedric model. The result 
indicates, that for k > only the ferromagnetically ordered low temperature 
phases remain stable, i.e. layered phases are thermodynamically suppressed at 
low but non zero temperature. We use Fade approximations to calculate crit- 
ical exponents of the magnetization, susceptibility and specific heat from the 
corresponding low temperature expansions in section 4. Finally we discuss the 
surface tension. There is some evidence, that a roughening transition occurs 
below the bulk critical temperature Tc. 



2 Ordered low temperature phases for k = 

If k vanishes, then only edges, where two plaquettes meet at a right angle cost 
energy. These edges are surrounded by an odd number of negative spins. The 
flipping of a whole (d — l)-dimensional spin layer therefore does not change the 
energy, i.e. apart from the global Z2 symmetry, the Hamiltonian shows an 
additional layer-flip symmetry at = 0. The ground states are the ferromag- 
netically ordered states and all configuration, which are connected with them 
by layer-flip operations. In the following we will show, that for each of these 
ground states there exists an ordered low temperature phase. These phases 
can be characterized by a non vanishing "spontaneous magnetization" , which 
for a finite system with spins is defined by 

M„ . (toL>. ,4) 

Here refers to a boundary condition, which singles out one ground state. 

denotes the number of spins, which are flipped compared to this ground 
state and accordingly the number of spins, which are unchanged. To show, 
that Mn does not vanish at sufficiently low temperatures, we will use a mod- 
ified Feierls argument PT|P^ . The idea is the following: Consider a finite 
system, where the spins at the boundary are fixed, such that one ground state 
is favoured (see figure |^). If we now swap a little "droplet" of spins inside 
the volume, the amount of energy we need will be essentially proportional to 
the number / of simple edges, which were established by swapping the spins. 
The simple edges form connected edge diagrams. Each overturned spin is sur- 
rounded with at least one such diagram. We can therefore estimate 
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Fig. 2. Starting from a ferromagnetic ground state, further ground states of the 
gonihedric model for A; = can be constructed by swapping (d — l)-dimensional 
spin layers. The figure above visualizes such a ground state. In order to carry out 
the Peier Is- argument, the boundary spins are fixed correspondingly. 

where ( 2^(^-1) ) ^ maximum number of spins an edge diagram with 

/ simple edges encloses and p{l) denotes the probability of occurrence of such 
a diagram. p[l) can be further estimated by 

p{l)<gil)e-^'\ (6) 



where g{l) is the entropy factor, i.e. the number of connected edge diagrams 
with I edges. We will show, that g{l) does not grow faster than exponentially. 
Thus (^) will be arbitrarily small for sufficiently large j3. The same type of 
argument was used in |^| to show the existence of a phase transition in the 



gonihedric model for A; > 0. The argument given there however contains a 
flaw, since the edge diagrams are not independent from each other, ii k > 0[^. 
We will now proceed with the details for the case k = 0. 

Consider a flnite system with spins in its interior and flx the spins at the 
boundary, so that they belong to a ground state, as indicated in flgure |^. 
For a conflguration cr of the spins inside the volume let C[(t] be the set of all 
(((i — 2)-dimensional) edges of the dual lattice, where two ((rf — l)-dimensional) 
plaquettes of the domain wall meet at a right angle. Each of these edges is 
surrounded by four spin, whose product is -1. Two edges in C[cr] are connected, 
if they have at least one ((ci — 3)-dimensional) vertex in common. C[(t] can be 
uniquely decomposed into connected parts: 

C[cT]=C[^U...UCl-, (7) 



^ For details see [||] 
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where C- denotes a connected component with li edges. The total energy E[a] 
can be written as 



E[cT]=J2Jh. (8) 



1=1 



Let g{l) be the number of connected edge diagrams with / edges. For each 
diagram we introduce a variable x', ^ = 1, ■ ■ ■ , 9(1), such that 



1 if the corresponding edge 

diagram occurs in cr (g) 

otherwise. 



Since each spin, which is flipped compared to the ground state, is surrounded 
with at least one connected edge diagram, we can estimate 



9(0 



Thus for the expectation value of the number of the overturned spins we 
obtain: 

/ / \^ 9(1) 

^ ? (2d(Z3T)j |:(xs>.. (11) 



If we number the spin configurations, which contain C' by «Ti,...,<Tfe, the 
expectation value {xDn can be written as: 



i^^N-^t^-''^"^^, (12) 

where 

^iv^E^-^'''"^- (13) 



For each (Tj there exists a unique configuration cr*, which does not contain Cl 
but leaves all other edge diagram unchanged, i.e. 

E[cr*] = E[a,]-Jl. (14) 
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The configurations cr*, . . . , cr^ are pairwise different. By restricting the parti- 
tion function to these configurations, we can estimate: 



k 



Zn>Y. e"^"""^' = e^^' E (15) 
i=i i=i 



Together with (|T^ this yields 



{x\)n <^-^''- (16) 



Hence the expectation value of the number of overturned spins is bounded by 



d 
d-2 



<*->"^Sb(^l P(Oe- (17) 



To proceed with the argument, we need an upper bound for g{l), the number 
of connected edge diagrams with / edges. The edges in those diagrams are 
connected via {{d — 3)-dimensional) vertices. What kind of edge configura- 
tions can occur at a vertex is determined by the possible configurations of the 
eight surrounding spins. As shown in figure |^, only four types of vertices can 
arise. With the following method all possible connected edge diagrams can be 

2- Vertex 3- Vertex 4- Vertex 5- Vertex 




Fig. 3. The possible edge configurations around a vertex are determined by the eigiit 
surrounding spins. Only four types of configurations can arise (modulo rotations) 
as shown above. 

constructed: 



1. Number all Vertices {{d — 3)-dimensional plaquettes) of the lattice. 

2. Choose one of the lattice points, e.g. r = (ri, . . . , r^), e Z and attach 
the edge k = {(ri, r2, A3, ... , A^) | r-j < Aj < rj -|- 1}. This is possible, since 
edges of all orientations occur in every possible diagram. 

3. Imagine that a connected subdiagram exists already. Consider all vertices 
in that diagram, where the configuration of the surrounding edges is not 
allowed (see figure H). Choose the one with the lowest number and attach 
further edges, such that an allowed configuration results. 

4. Repeat constructions step 3, until all I edges are attached. 
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The maximum number of outcomes of the procedure given above is an upper 
bound for g{l). For construction step 2 there are N possibihties. Each time we 
perform constructions step 3, at least one edge is added. The number of choices 
to complete the vertex depends on the edge configuration already present. In 
figure ^, we list all possibilities which can arise together with the number of 
choices n{k) to attach k edges. Let n^^^{k) be the maximum number of choices 
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Fig. 4. The figure shows all possibilities of incomplete vertices (modulo rotations). 
n{k) is the number of choices to complete the vertex by adding k edges. 



to attach k edges. The function 

fr{x):=Nx(^j2n^^Mx'''^ , (18) 



can now be interpreted as a generating function, which counts all connected 
edge diagrams with / edges at least once, which can be constructed by repeat- 
ing constructions step 3 r-times, i.e. the coefficient of in the series expansion 
of fr{x) is an upper bound for the number of those diagrams. Therefore an 
upper bound for g{l) is provided by 



9{l)< 




(19) 



x=0 



The values n^^^{k) can be read off from figure |^. We find: 



1 



dx 



Nx 



1 - (2a; + 4a;2 + 2x3 + x^ + x^) 



x=0 



G{x) 



I z=0 ' 



(20) 
(21) 
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G{x) can be written as: 



G{x) 



Nx 

1 - (2a; + 4a;2 + 2x^ + x^ + x^ 



(22) 



Putting this expression in (pO|), we obtain 



\ XkJ \XkJ 



(23) 



The increase of g{l) is dominated by the smallest \xk\ ='■ x^^^. Hence we can 
further estimate: 



g{l) <N \ < Aro.67(3.39)^ 

\A;=1 J V Fm| / 



(24) 



Substituting this result in (|3), we find an upper bound for the density of 
overturned spins: 



N 



<E 



/ \2d{d-l)^ 



0.67(3.39)' e 



(25) 



This inequality remains valid in the thermodynamical limit. The sum on the 
right hand side approaches zero, if /? tends to infinity. Therefore, for (3 > 13c it 
will be smaller than which implies a non-zero spontaneous magnetization. 
This completes the proof. Note that the estimate above is only valid for d > 
3. Indeed in two dimensions, the model for A; = is equivalent to the 1- 



dimensional ordinary Ising model |T2|. Thus no transition occurs in this case. 



The upper bound for (iV_)^ given above can be used to calculate an upper 
bound for Pc- We find 



PcJ < 1.43 for d = 3, 
(3cJ <12?, for d = A. 

The upper bound for f3cJ in 3 dimensions is consistent with values obtained 
from simulations p6| , p7| {j3cJ ~ 0.505), cluster variation-Pade approximations 
H {(5cJ ^ 0.550) and mean field calculations |2§ {p^J ~ 0.325). 



If A; > a similar proof can not carried through, since in this case the connected 
edge diagrams are not independent from each other, i.e. given a configuration 
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cr which contains a certain edge diagram C|, the corresponding configuration 
cr* does not necessarily fulfil the condition 

E[a*] = E[a]~E[C% (26) 



where -E[C'] denotes the energy contribution of the connected edge diagram. 
In going over from cr to cr*, the energy can even increase by an amount ~ t^. 



3 Low temperature expansion of the free energy 



If A; > then the Hamiltonian (|^) is no longer invariant under layer-fiip oper- 
ations. Only the global Z2 spin-fiip symmetry remains. However fiat domain 
walls still do not cost energy. Thus apart from the ferromagnetically ordered 
ground states, there are layered ground states containing fiat domain walls 
which do not intersect. As was already mentioned in [^, the ground state 
contributions Z^^^ of the total partition function 

Z = Y, Z^3^ (27) 



are not degenerate in contrast to the A; = case. Here Z^^'' contains the ground 
state g and all configurations, which differ from ghy a. "small" number of spins. 
The lowest excitation is realized by a configuration, where one spin is fiipped 
compared to the ground state. Since this configuration contains 12 simple 
edges (elementary cube), its energy is 12 J, independent of the underlying 
ground state. The next higher excitation is obtained by swapping an additional 
nearest neighbour spin. The energy of this configuration however depends on 
the ground state: If there is no fiat domain wall present between those two 
spins, the energy is 16J. Otherwise four additional double edges contribute, 
i.e. the energy amounts to 16J + 16Jk. Thus to lowest order, layered low 
temperature phases seem to be energetically disfavoured |2^. In this section 
we will quantify this effect. For this we perform a low temperature expansion 
of the free energy for the three dimensional case to estimate the magnitude 
of each ground state contribution Z^^\ This analysis indicates, that at low 
but non-zero temperature the occurrence of layered low temperature phases 
is indeed thermodynamically suppressed. 

Consider a finite system of volume V = L^. By N we denote the set of the 
N = Lr" spins inside V . Each excitation of a given ground state g is defined by 
the subset I C N of spins, which are swapped compared to g. The excitation 
energy E'^^^il) depends on g. We define the ground state contributions Z^-* 
by: 
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ICN 

where 

,i9)(l):=e-(^^'''('\ (29) 

Ground states g which can be mapped onto each other by a rotation or a 
global spin flip yield the same Zj^\ We therefore write: 

Z^ = 2(4°)+3'e"'4')), (30) 

9=1 



where g refers to ground states modulo rotations and global spin-flip and g = 
denotes the ferromagnetically ordered ground state. The logarithm of Z^'' can 
be written as: 



ln(4))= 5:c(^)(I), (31) 

ICN 



where the cluster contributions c'^^^(I) are defined by 

c(3)(I) = ^(_i)(lil-IJ|) ln( J2 z^'\K)). (32) 



JCI KCJ 



As can easily be shown, these cluster contributions vanish, if the cluster I is 
not "connected" . Here we call two spins connected, if their elementary cubes 
have at least one edge in common. In general two spin clusters Ii,l2, which 
can be mapped onto each other by a translation will contribute differently, 
unless the translation is parallel to the ground state planes. To take care 
of this restricted translation invariance we introduce structural coefficients, 
which characterize the ground states such that the right hand side of (|HT|) 
can be expanded systematically: For each ground state g there exists a unique 
(L — l)-tupel p{g) = {pi, . . . ,pl_i), pi = 0, 1, which characterizes the sequence 
of ground state planes say in positive direction. Pi = 1 means, that a ground 
state plane is present (see figure H). For g = (gi, . . . , g„), gj = 0, 1, <= L — 1 
and 

T := {k eN \ qi = pk,. . . ,qn= Pk+n-i}- (33) 



we define: 



Qq := Number of elements in T. (34) 
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L = 10, = (0,0, 1,0, 1,0, 0,0,1) 

Fig. 5. The picture above displays ground state layers within a finite volume of 
linear length L = 10. The sequence of this layers in positive direction (from left to 
right) is described by a (L — l)-tupel p. 

Thus aq{g) is the number of translations, such that the sequence q matches 
with a subsequence of p{g). If q does not have any elements, we set a{g) = 
L. The structural coefficients are not independent from each other, e.g. the 
following equations hold: 

a{g) = L 
ao{g) + ai{g) = L - 1 

(35) 

dooig) + aoi{g) + aioig) + «ii(fl') = L-2 



Using the structural coefficients, we can now write: 

lim hm (a(^)/ + ao(y)/o + ai(^)/i + aoo(^)/oo + ---)-(36) 

L—too iV L— »oo J_/ ^ ' 



In this expression / denotes the sum of all cluster contributions of widths 1, 
i.e. all contributions from clusters, whose linear extension perpendicular to 
the ground state planes is 1. /o denotes the sum of all cluster contributions of 
widths 2, which are not cut by a ground state plane, /i denotes the sum of all 
cluster contributions of widths 2, which are cut by a ground state plane and 
so on. A systematic low temperature expansion of the cluster functions /g can 
now be generated by analysing the corresponding cluster contributions. We 
obtain 
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JCI KCJ 



(37) 



oo f -i\ (r— 1) r 

= E^^ E n-^^HK,) 

r=l ' KiU...UK,.=I j=l 

K,CI, K,^0 

= ^^'\l)-l E ^^^HKi);^(^)(K2) + ... (38) 

^ KiUK2=I 
Ki,K2CI 

Ki,K2^0 

The Boltzmann factors z^^^ (I) are of the form: 



^(^)(I) = x"^l/"*, (39) 

where x = e~^'' , n2 = Number of simple edges 
y = e~^'^^-^ , 77-4 = Number of double edges. 

From here one can see, that a cluster I contributes to an order in x, which is 
higher or equal to 4(6i(I) + 62(1) + ^3(1)), if hij) denote the side lengths of 
the smallest box, which contains I. Thus to generate an expansion in x up to 
order n, only clusters I with 4(6i(I) +62(1) +^3(1)) < have to be considered. 

With the method described above, we calculated the cluster functions fg up 
to order x^^. The numbers of clusters which contribute to a specific order are 
listed in table |l|. To see, whether or not layered low temperature phases arise, 
we write according to (pO|): 

7 2^-1-1 7(9) 

^ = 1 + 3 E % (40) 



97W ^ 7, , 



Using (|36|) we arrive at: 



2^-1-1 



lim — ^ = 1 + 3 lim ^ expi L^(ai{g){fi - /o) + am{g){fm - /oo) 

+aio(^7)(/io - /oo) + a^MUii " /oo) + . . .)). (41) 

Since the functions (/i — /o),(/oi ~/oo)v . . turn out to be negative in the range 
< X < 1 and vanish only at a; = and x = 1 for all positive fc, the right hand 
side of (Hll) tends to 1 in the thermodynamical limit. It follows that layered 
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"\T 1 r * 1 1 

M umber oi simple edges 


"\T 1 r 1 J 

M umber oi clusters 


12 


1 


14 





16 


3 


18 





20 


6 


22 


18 


24 


10 


26 


96 


28 


105 


30 


372 


32 


789 


34 


1806 


36 


4881 


38 


10134 



Table 1 

Numbers of connected spin clusters with a given number of simple edges. 

phases are suppressed at low but non-zero temperature. In figure ^ the function 
(/i ~ fo)/^^^ is plotted for k = 1. Higher order functions (/oi — foo)/x'^'^,. . . 
show the same qualitative behaviour. The total free energy F can therefore 




Fig. 6. (/i — fo)/x^^ is plotted as a function of x = e ^ for k = I. 

be calculated as a low temperature expansion, where only excitations of the 
ferromagnetic ground states are considered, i.e. 

-PF = ln(2) + ln(Z(0)) = ln(2) + f + f, + f,, + f + . (42) 
To order x^^ this yields: 
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-f3F = ln(2) + + 3x^^ + 6x^° + {6y + 12)x'^^ + ^x^^ + (48 + 48y)x'^^ 

+ (15 + 6y^)x^^ + (136 + 204y + 2Ay^ + 8?/^)x^° + {228y + QOy^ 

129 1975 
+— + (456y + 216y2 + 78y^ + 456)x^^ + (1824?/ + — - 

+1032?/2 + 56?/3 + 2Ay^ + 2?/^)x3'^ + (1110?/ + 1224?/^ + 600?/^ 
+108?/^ + 540 + 2Ay^)x^^ + ... 

where x = e"^*^ and y = q-P^^J , With the formalism described above, it is 
also possible to derive a low temperature expansion of the surface tension cr, 
which is defined as: 



a:=lim — In^ (43) 

= (/o - /i) + (/oo - /oi) + (/oo - /lo) + (/ooo - /ool) + • • • (44) 

In this expression Z^'* denotes the contribution of the total partition function 
from a ground state which contains exactly one plane domain wall. Clearly a 
vanishes for T = 0. However for finite T a finite surface tension is generated 
by thermal fluctuation. We flnd the following series expansion: 



a = (1 - y^)x^^ + (-2?/^ - 2y^ + 4)x=^° + (4?/ + 4 - 8y^)x^'^ + (-4?/^ (45) 
+10 - 3?/^ - 3y^)x^^ + (16 + 48?/ - 8y^ - 2Ay^ - 2Ay^ - 8y^)x^^ 
+ (-10?/*^ + 72 - 32?/3 - 16?/^ - 6?/^° - 18?/2 + 16?/^ - Qy^)x^^ 
+ (32 + 264?/ + 16?/2 - 40?/^ - 48?/^ - 40?/^ - 120?/^ - 16?/* - 48?/^)x^° 

65 

+ (240?/ + 4?/2 - 200?/^ - 27 y^ - 4?/^ - 96?/^ - —y^ - 32?/^° - 40?/^ 

-13?/^2 + ^)x=^2 + (828?/ + 104?/2 + 48?/^ - 304?/^ - 160'i/^ - 540?/^ 

-136?/* - 24?/^° - 296?/^ - 96?/^ - 4?/" + 580)x^^ + (2208?/ + 1114?/^ 
-792?/^ - 718?/^ - 344?/^ - 928?/^ - 146?/* - 88?/^° - 424?/^ - 90?/^^ 
-64/ - 16?/^^ - 26?/^^ + 314)x^^ + (3172?/ + 568?/^ + 632?/^ - 1056?/^ 
-1272?/^ - 1268?/^ - 712?/* - 352?/^° - 1296?/^ - 64?/^^ - 748?/^ 
-232?/^^ - 16?/^^ + 2644)x^* + . . . 



4 Critical exponents 

If the self-avoiding coupling k is sufficiently large {k > 0.3), the authors of 



SB|,|27[| find a second order phase transition, whereas at /c = the transition 
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is of first order. Tliese results were obtained by simulations. Using the result 
above, i.e. that only the ferromagnetically ordered phases contribute to the 
free energy, we were able to derive low temperature expansions for the mag- 
netization, susceptibility and specific heat up to order for fixed values of 
k {k = |, n = 0, . . . , 12). A Pade analysis yields critical exponents, which are 
in agreement with previous results. 



In the following we describe the method used to derive the series expansions 
35| , |36[| . The free energy can be derived from the ferromagnetic part Z^^^ of 



the partition function. According to (31) we write: 



HzS^) = E c^^Hi)- (46) 

ICN 



To each cluster I of overturned spins corresponds a minimal box q with side 
lengths qi,q2,(l3, which contains I. We define 

C{q) := Sum of all cluster contributions c^°)(I) modulo transla- (47) 
tions, such that q is the minimal box of I. 



If Z{Q) denotes the partition function for the finite box Q, evaluated for a 
ferromagnetic ground state, the following equation holds due to translation 
invariance: 

ln(Z(g)) = J2(Q,-q, + 1)(Q2 - q2 + 1)(Q3 - ^3 + l)C(g). (48) 

qCQ 



Therefore the free energy can be written as: 

-PF = J2Ciq)= E C(g) + 0(x"+2). (49) 

Q q : 4:{qi+q2+q3)<=n 



Inverting equation ( P5| ) and putting the resulting expression for C{q) in 
we finally arrive at: 



(3F= E (-l)'-«i-'32-Q3 in(Z(g)) (50) 



+0{x 

where = 0, ifn< Oorn> 5. Thus in order to generate a low tem- 
perature expansion for F, we just need to calculate partition functions for 
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small volumes, which can be done very efficiently by recursive methods^ ||37|| . 
Expanding 



Z{Q) 



. ^ ^n2{I)+4fen4(I) ^n(I) 
ICN 



Z 



Kf\x) + - 1) + -lf + .. 



(Q) 



(51) 
(52) 



where n{I) denotes the number of spins in I, series expansions for the magne- 
tization M{x), susceptibility S{x) and specific heat C{x) can be derived from 
the polynomials Aq'^''(x), A^'^''(x),A2'^''(a;) in the following way: 



M{x) 
S{x) 
C{x) 



1-2 



d{-PF) 



dz 



z=l 



^ d f d{-pF) 

dz ^ dz 
x—lx 



dx 



dx 



z=l 



2 = 1 



The resulting expansions for fixed values of A; = ^ 



4 ' 



n 



0, 



(53) 
(54) 
(55) 

, 12 can be found 

in appendix A. To calculate critical exponents, we analysed [n/ (n — 2)], [n/{n — 
1)], [n/n], [njin + 1)] and [n/(n + 2)] unbiased Dlog Pade approximants of 
the corresponding series expansions [Q. The results are shown in table ^, ^ 
and |. The errors roughly estimate the fluctuations within the Pade tables. 
Defective approximants have been excluded from the analysis. In the 
critical x was determined as Xc ~ 0.64. A cluster variation-Pade calculation 
also confirms this result The authors determine the critical exponent of 
the magnetization as /? = 0.062 ± 0.003 at a critical temperature Xc ~ 0.65. 
The exponent 7 was also calculated in |^ as 7 ~ 1.6. An estimate of 7 ~ 1.4 
was given in |30|. Our values are of the same order of magnitude. However, 



the fluctuations within the Pade tables are relatively large. Unfortunately 
different extrapolation techniques like inhomogeneous Pade approximants or 
differential approximants did not yield better results. The exponent a of the 
specific heat was determined in I^Sj as a ~ 0.7. Our values agree roughly 
with this result. Since the errors are quite large also in this case, it is unclear 
whether the exponents are fc-dependent or not. Furthermore we analysed the 
series expansion of the surface tension a{x,y) as given in (^). The surface 
tension should scale at the critical point and vanish for higher temperatures. 
Surprisingly, the Pade analysis of the series expansion gives negative values 



A detailed description of the used algorithm can be found in 
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k 


X(. 


P 


0.50 


0.636 ± 0.002 


0.037 ±0.002 


0.75 


0.642 ± 0.003 


0.045 ± 0.004 


1.00 


0.641 ± 0.001 


0.040 ± 0.002 


1.25 


0.644 ±0.001 


0.045 ± 0.002 


1.50 


0.646 ±0.001 


0.047 ±0.002 


1.75 


0.647 ±0.002 


0.049 ± 0.003 


2.00 


0.647 ±0.002 


0.050 ± 0.004 


2.25 


0.647 ±0.002 


0.048 ± 0.007 


2.50 


0.648 ± 0.003 


0.050 ±0.006 


2.75 


0.648 ± 0.002 


0.051 ±0.006 


3.00 


0.649 ± 0.002 


0.052 ± 0.004 



Table 2 

Critical temperature Xc = e~^'"^ and magnetization exponent /?. 



k 




7 


0.50 


0.642 ± 0.002 


1.7±0.1 


0.75 


0.647 ±0.008 


1.8 ±0.3 


1.00 


0.64 ±0.01 


1.7 ±0.2 


1.25 


0.65 ±0.01 


1.7 ±0.3 


1.50 


0.65 ±0.01 


1.7±0.3 


1.75 


0.65 ±0.01 


1.6 ±0.3 


2.00 


0.65 ±0.01 


1.7 ±0.3 


2.25 


0.64 ± 0.02 


1.7 ±0.3 


2.50 


0.65 ±0.02 


1.7 ±0.7 


2.75 


0.64 ±0.02 


1.5 ±0.8 


3.00 


0.65 ± 0.02 


1.3 ±0.6 



Table 3 

Critical temperature Xc = e~l^'"^ and susceptibility exponent 7. 
for the corresponding exponent /i. Moreover, a{x,y) can be written as: 

cr{x,y) = {1 -y)P{x,y), (56) 
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k 




a 


0.50 


0.647 ±0.002 


0.83 ±0.03 


0.75 


0.646 ±0.004 


0.67 ±0.08 


1.00 


0.65 ±0.01 


0.62 ±0.03 


1.25 


0.649 ±0.001 


0.57 ±0.01 


1.50 


0.649 ±0.001 


0.53 ±0.01 


1.75 


0.648 ±0.002 


0.49 ±0.02 


2.00 


0.648 ±0.001 


0.45 ±0.02 


2.25 


0.64 ±0.03 


0.5 ±0.1 


2.50 


0.647 ±0.003 


0.41 ±0.06 


2.75 


0.64 ±0.03 


0.4 ±0.2 


3.00 


0.63 ±0.03 


0.3 ±0.2 



Table 4 

Critical temperature Xc = e~^'"^ and specific heat exponent a. 



where all coefficients in the expansion of P{x,y) are positive. This indicates 
a diverging surface tension, which is of course not the expected behaviour. 
However a similar situation arises in the ordinary 3-dimensional Ising model. 
The low temperature expansion of the surface tension [R^] reads in this case: 



242 

a{x) = -ln(x) - 2 - 2 x*^ - 10 - 16 x^^ - 150 

3 



-734 x^^ - 



4334 



■ X 
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X 



-2/3 J 



)• 



Here the expansion coefficients all have the same sign and a Fade analysis also 
yields negative critical exponents. Nevertheless, simulations indicate, that the 



surface tension scales at the critical point [|^. The corresponding exponent 
is positive and consistent with Widom-scaling. The reason for this strange 
behaviour lies in the existence of a roughening transition at a temperature Tr 
below Tc EO]. Renormalization group calculations indicate, that this transition 



is of Kosterlitz-Thouless type ||5l[] with an essential singularity at Tr. We take 
this as a strong hint, that a similar transition occurs in the surface tension of 
the gonihedric Ising model. This point needs however further clarification. 
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5 Summary 



The gonihedric Ising model was introduced as a lattice discretization of the so 
called gonihedric string. The gonihedric string itself is a model for triangulated 
random surfaces. Whether or not both theories are equivalent in the sense 
that they are different discretizations of the same continuum random surface 
theory is not obvious. A necessary condition for the existence of a continuum 
limit is the occurrence of a second order phase transition. The existence of 
ordered low temperature phases often signals such a phase transition, if at 
high temperatures a disordered phase is expected. 

In this paper we studied some low temperature properties of the gonihedric 
Ising model. U k = 0, the model shows apart from the global Z2 symmetry an 
additional layer-flip symmetry. With a suitable extension of Peierls contour 
method we proved, that the layer-flip symmetry is spontaneously broken at 
low temperatures for d > 3. We found infinitely many ordered low temper- 
ature phases, one for each ground state. However, simulations indicate, that 



the corresponding phase transition is of first order |27|. Thus in this case a 



continuum limit does not exist. For k > 0.3 the authors of [0,EB,E3,Eni find 



a second order phase transition. It was suggested in I^Q], that in contrast to 
the case k = only the ferromagnetic phases remain stable, while layered 
phases are thermodynamically suppressed. The low temperature expansion of 
the free energy in section 2 supports this picture. For k > 0, a finite surface 
tension is generated if T > 0, i.e. the occurrence of layers is energetically sup- 
pressed by a factor e~^^ . On the other hand the ground state entropy grows 
only like 2^ with the linear size L of the system and can therefore not com- 
pensate the energy factor. Using this result, we calculated low temperature 
expansions of the magnetization, susceptibility and specific heat, where only 
excitations of the ferromagnetic ground states were considered. The critical 
exponents, derived from Fade approximants are in agreement with previous 
results. Unfortunately the fluctuations within the Fade tables are quite large 
(in particular for 7 and a). It is therefore not possible to determine, whether or 
not the model shows non-universal behaviour. Finally we analysed the surface 
tension. A Fade- approximation resulted in a negative critical exponent. This 
probably indicates the existence of a roughening transition as in the ordinary 
3-dimensional Ising model. 



A Series expansions 



This appendix contains the coefficients of the series expansions, used in section 
4 to calculate critical exponents. The corresponding Fade tables can be found 
in pi. 
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A . 1 Magnetization 



k 





0.25 


0.5 


0.75 


1 


1.25 


1.5 




1 


1 


1 


1 


1 


1 


1 




-2 


-2 


-2 


-2 


-2 


-2 


-2 




-12 


-12 


-12 


-12 


-12 


-12 


-12 




-42 


-42 


-42 


-42 


-42 


-42 


-42 




-96 


-72 


-72 


-72 


-72 


-72 


-72 


x23 





-24 

















x24 


-74 


-74 


-98 


-74 


-74 


-74 


-74 


^25 











-24 











x26 


-720 


-432 


-432 


-432 


-456 


-432 


-432 


x27 





-288 











-24 





x28 


-564 


-516 


-804 


-516 


-516 


-516 


-540 


^29 











-288 











x^o 


-3680 


-1736 


-1688 


-1688 


-1976 


-1688 


-1688 


^31 





-1752 











-288 





^32 


-5316 


-2856 


-4464 


-2664 


-2664 


-2664 


-2952 







-1920 





-1752 













-18048 


-9468 


-10752 


-8736 


-10440 


-8688 


-8688 


^35 





-6456 





-1872 





-1752 





^36 


-46046 


-17582 


-22610 


-15518 


-17246 


-15326 


-17078 







-20184 





-6456 





-1872 





^38 


-84696 


-41472 


-52992 


-31980 


-37848 


-31248 


-33072 







-298-iO 





-19581 





-G456 







-337326 


-118098 


-141546 


-103698 


-121758 


-101634 


-107946 


x^i 





-128760 





-29184 





-19536 





x42 


-518832 


-236876 


-281280 


-152344 


-173560 


-142852 


-161800 


^43 





-174816 





-122784 





-29184 







-2016750 


-619662 


-739530 


-496470 


-612462 


-482310 


-510018 


^45 





-731672 





-159560 





-122184 






k 


1.75 


2 


2.25 


2.5 


2.75 


3 


X^ 


1 


1 


1 


1 


1 


1 


Xl2 


-2 


-2 


-2 


-2 


-2 


-2 




-12 


-12 


-12 


-12 


-12 


-12 


x20 


-42 


-42 


-42 


-42 


-42 


-42 


x22 


-72 


-72 


-72 


-72 


-72 


-72 


x24 


-74 


-74 


-74 


-74 


-74 


-74 


x26 


-432 


-432 


-432 


-432 


-432 


-432 
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^28 


-516 


-516 


-516 


-516 


-516 


-516 


^29 


-24 

















^30 


-1688 


-1712 


-1688 


-1688 


-1688 


-1688 


^31 








-24 











^32 


-2664 


-2664 


-2664 


-2688 


-2664 


-2664 


^33 


-288 











-24 







-8688 


-8976 


-8688 


-8688 


-8688 


-8712 


^35 








-288 











^36 


-15326 


-15326 


-15326 


-15614 


-15326 


-15326 


^37 


-1752 











-288 





^38 


-31200 


-32952 


-31200 


-31200 


-31200 


-31488 


^39 


-1872 





-1752 











^40 


-101442 


-103314 


-101442 


-103194 


-101442 


-101442 




-6456 





-1872 





-1752 







-142120 


-148528 


-142072 


-143944 


-142072 


-143824 




-19536 





-6456 





-1872 







-18021(3 


-199(338 


-180051 


-18(3510 


-180051 


-181926 




-29184 





-19536 





-6456 






A. 2 Susceptibility 



k 





0.25 


0.5 


0.75 


1 


1.25 


1.5 




1 


1 


1 


1 


1 


1 


1 




12 


12 


12 


12 


12 


12 


12 




75 


75 


75 


75 


75 


75 


75 




132 


108 


108 


108 


108 


108 


108 


x23 





24 

















x24 


290 


290 


314 


290 


290 


290 


290 


x25 











24 











x26 


1416 


984 


984 


984 


1008 


984 


984 


x27 





432 











24 





x28 


2274 


2178 


2610 


2178 


2178 


2178 


2202 


x29 











432 











^30 


9664 


5500 


5404 


5404 


5836 


5404 


5404 







3804 











132 





X32 


20064 


14742 


18258 


14358 


14358 


14358 


14790 


^33 





3936 





3804 












22 



^34 


65916 


39138 


41544 


37392 


41100 


37296 


37296 







21324 





3864 





3804 





^36 


170779 


93603 


110913 


88059 


91635 


87675 


91479 


^37 





54192 





21324 





3864 





^38 


420492 


233664 


264696 


207786 


227652 


206040 


209808 


^39 





138896 





52896 





21324 







1400784 


660426 


770154 


610026 


658764 


604482 


625518 







445740 





136968 





52824 







2963424 


1526998 


1701136 


1244228 


1359476 


1218350 


1269716 


^43 





983064 





428340 





136968 







9979599 


3989439 


4711941 


3507195 


3906159 


3457395 


3590277 







3297892 





935276 





427044 






k 


1.75 


2 


2.25 


2.5 


2.75 


3 




1 


1 


1 


1 


1 


1 


^16 


12 


12 


12 


12 


12 


12 




75 


75 


75 


75 


75 


75 




108 


108 


108 


108 


108 


108 


OA 


290 


290 


290 


290 


290 


290 




984 


984 


984 


984 


984 


984 


^28 


2178 


2178 


2178 


2178 


2178 


2178 




24 

















^30 


5404 


5428 


5404 


5404 


5404 


5404 


X^^ 








24 











^32 


14358 


14358 


14358 


14382 


14358 


14358 




432 











24 





^34 


37296 


37728 


37296 


37296 


37296 


37320 


^35 








432 











^36 


87675 


87675 


87675 


88107 


87675 


87675 


^37 


3804 











432 





^38 


205944 


209748 


205944 


205944 


205944 


206376 


^39 


3864 





3804 













601098 


607962 


00 1098 


607902 


601098 


001098 




21324 





3864 





3804 







1216604 


1237832 


1216508 


1220372 


1216508 


1220312 




52824 





21324 





3864 







3451851 


3504387 


3451467 


3472791 


3451467 


3455331 




136968 





52824 





21324 
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A. 3 Specific heat 



k 





0.25 


0.5 


0.75 


1 


1.25 


1.5 




144 


144 


144 


144 


144 


144 


144 




768 


768 


768 


768 


768 


768 


768 




2400 


2400 


2400 


2400 


2400 


2400 


2400 




8712 


5808 


5808 


5808 


5808 


5808 


5808 







3174 



















288 


288 


3744 


288 


288 


288 


288 













3750 













64896 


32448 


32448 


32448 


36504 


32448 


32448 







34992 











4374 





^28 


16464 


11760 


49392 


11760 


11760 


11760 


16464 


^29 











40368 













334800 


127800 


122400 


122400 


165600 


122400 


122400 







196044 











46128 





^32 


391680 


90624 


281088 


66048 


66048 


66048 


115200 


^33 





257004 





222156 











^34 


1394136 


631176 


818448 


534072 


762960 


527136 


527136 


^35 





558G00 





279300 





219900 





^36 


4660848 


1133136 


1571184 


884304 


1156464 


853200 


1117584 







2603838 





624264 





312132 





^38 


5207064 


2269968 


3725520 


909720 


1472880 


788424 


1108992 


^39 





1773486 





2786472 





693576 







36662400 


9576000 


11131200 


7924800 


10641600 


7617600 


8318400 







15956052 





1865910 





3066144 







36694728 


15734880 


20095488 


3972528 


4505256 


2310840 


5411952 







10472736 





16585530 





2052390 







213025824 


52614672 


58428480 


36509088 


53352288 


34557600 


36462624 







80842050 





8885700 





18022500 






k 


1.75 


2 


2.25 


2.5 


2.75 


3 




144 


144 


144 


144 


144 


144 




768 


768 


768 


768 


768 


768 




2400 


2400 


2400 


2400 


2400 


2400 


99 

X 


5808 


5808 


5808 


5808 


5808 


5808 




288 


288 


288 


288 


288 


288 



24 



^26 


32448 


32448 


32448 


32448 


32448 


32448 


^28 


11760 


11760 


11760 


11760 


11760 


11760 


^29 


5046 

















^30 


122400 


127800 


122400 


122400 


122400 


122400 


^31 








5766 











^32 


66048 


66048 


66048 


72192 


66048 


66048 


^33 


52272 











6534 







527136 


582624 


527136 


527136 


527136 


534072 










58800 











^36 


853200 


853200 


853200 


915408 


853200 


853200 




279276 











65712 





^38 


779760 


1074336 


779760 


779760 


779760 


849072 


X^^ 


346788 





310284 













7579200 


7944000 


7579200 


7905600 


7579200 


7579200 


X« 


766536 





383268 





342924 







2162664 


2956464 


2152080 


2554272 


2152080 


2511936 




3372570 





813111 





-121572 







34185888 


37682304 


34139424 


35022240 


34139424 


34580832 




2247750 





3693600 





923400 
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